Theory of Orbital Nematicity in Underdoped Iron Arsenides 
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Recent finding of an unusual in-plane resistivity anisotropy in the underdoped 122-family at high 
temperature (T) suggests an orbital nematic (ON) order, posing a challenge to extant theories. 
The sign of the anisotropy contradicts expectations from weakly correlated as well as pure spin-only 
nematic views. Here, we show how such an ON order with accompanying structural distortion arises 
from residual, intersite and inter-orbital two-body interactions in an incoherent "bad metal" close to 
Mottness. Enhancement of orbital-selective incoherence is shown to be necessary for understanding 
transport anisotropy. Our results suggest that ON order, with subsequent antiferromagnetic order 
might be the primary competitor to superconductivity in Fe arsenides. 

PACS numbers: 74.70.Xa, 74.25. F- 74.20.Mn 



The "Iron Age" of high temperature superconductiv- 
ity (HTSC) appears to reinforce the paradigm shift set in 
motion by HTSC in cuprates more than twenty years ago. 
Examining analogies and differences between cuprates 
and Fe-pnictides (FePn) can shed more light on both. 
While parent cuprates are antiferromagnetic (AF) Mott 
insulators, parent FePn are multi-band AF metals, with 
the Fermi surface (FS) constituted of several small pock- 
ets. Notwithstanding this difference, the normal metallic 
state (by chemical doping) in both defies even a quali- 
tative description within traditional Landau Fermi Liq- 
uid (LFL) concepts. FePn thus appear to be anomalous 
metals giving way to unconventional superconductivity 
at high temperature (T) [l[. 

AF in FePn is universally preceded by, or coincides 
with, a tetragonal-to-orthorhombic (T-O) structural dis- 
tortion. Recent work on 122- 2] and 1111-based FePn 
reveals signs of enhanced nematic (N) susceptibility (xn) 
and electronic crystallinity upon underdoping, deepening 
the analogy with cuprates. Transport data puts strong 
constraints on theories purporting to describe the N- 
phase. In particular (i) inspite of the fact that the 
N order leads to b < a (here, a, b are planar unit cell 
distances), and (ii) the q = (tt, 0) AF implies antifcrro- 
magnetically aligned spins along a and ferromagnetically 
aligned spins along b, the resistivity Pb(T) along b exceeds 
p a {T) along a, an inexplicable result in a weak-coupling 
LFL [3] or pure spin- fluctuation [j| scenarios. In addi- 
tion, Pb(T) shows insulator-like behavior above T Cl even 
as p a {T) remains metallic. These observations call for a 
more careful investigation into microscopic sources of ne- 
maticity, especially in view of the resistivity anisotropy 
peaking (maximal nematic response, xn) around a dop- 
ing where SC T c (x) is highest. Thus, the incoherent "nor- 
mal" state gives way to either (a) a nematic phase (at 
least in 122-FePn), with a T-0 distortion and subsequent 
SDW order, upon underdoping, or (b) to an unconven- 
tional SC with a dome in T c versus x, or (c) to a renor- 
malized Fermi liquid metal at low T upon overdoping. 



These findings in 122-FePn conflict with pure spin fluc- 
tuation models ((i) and (ii) above), and, in addition to 
(i),(ii), in view of the insulator-like Pb(T), weakly cor- 
related models as well. Orbital-based models do iden- 
tify the T-0 distortion with multi-orbital (xz,yz) cor- 
relations @, but nematicity has hitherto not been con- 
sidered there. Finally, how such a coupled nematic-plus 
T-0 transition could arise as an instability of the "high"- 
T "strange" metal [H0| remains open (the T-0 transition 
generically precedes AF-spin density wave (SDW)). Here, 
we address these issues in detail using LDA+DMFT. Our 
philosophy is to identify the germs of the emergent or- 
der in the residual (two-particle) interactions present in 
the high-T incoherent fluid. Thereupon, extending our 
earlier LDA+DMFT calculations to a postulated orbital- 
nematic (ON) state, we show how the unusual features 
above find natural explication within a theoretical pic- 
ture where an orbital nematic order arises as a direct 
instability of a sizably correlated bad metal. 

In an incoherent metal, inter-band one-electron mix- 
ing is irrelevant in a renormalization-group (RG) sense, 
akin to the situation that occurs in weakly coupled 
D = 1 Luttinger liquids @. This inhibits a smooth 
crossover to a correlated FL, simultaneously favoring or- 
dered states arising directly from the incoherent metal. 
The underlying reason is that when one-electron mix- 
ing is irrelevant, the corresponding second-order pro- 
cess, corresponding to two-particle intersite and inter- 
orbital hoppings, turns out to be the most relevant. 
In our earlier LDA+DMFT for the "normal" state, an 
incoherent metal arises from an Anderson orthogonal- 
ity catastrophe-like effect due to orbital-selective (OS) 
Mott physics. The general two-body interaction, ob- 
tained at second order from a (now RG-irrelevant) one- 
electron inter-band term, t ab J2<i,j>,a( c iaa c jb<j + h.c), 

is = 5E a ,6,fe,fe'K6(fc,fc04,fe,t C 6,-fc,4. c 6-fc',4- c a,fc',t) 

where a, b = xz,yz and the scattering vertex is 

V a b(k,k',uj) = g 2 Xab(k - k',uj) with Xab(k,uj) being the 

inter-orbital susceptibility. The static, nearest- and next- 
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nearest neighbor parts of V a b(k, k') are V^ b \k,k',0) 



0(200 - 250) meV and V c 
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U' + Jh — " u "l LL ^ y " liuL "ab — U' + Jh — 

0(120 — 160) meV, consistent with magnetic excitation 
energy scales from INS data. These residual interactions 
contain coupled, inter-orbital charge and spin fluctua- 
tions in particle-hole (p-h) and cooper (p-p) channels, 
immediately suggesting the possibility of both, mutually 
competing instabilities Q from the "high" T incoherent 
metal. In earlier work, we showed how an unconventional 
SC with the right gap symmetry indeed arises from the 
cooper instability above 10 . Here, we study how a p-h 
decoupling of H^ 2 > above in the situation specific to FePn 
leads to the ON instability. We decouple in the p-h 
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12a,b,k,a(^b ( fc > k ') C i,k,a C b,k',a + 



sector as H 

b), with A^(fc,fc') = |MM^)(4 ife _, CT c 6>fc> ). Con- 
sistent with lattice symmetries, we write V a b(k, k') = 
^2iV^ b T}i(k)r}i(k'), with i]i(k) being the irreducible rep- 
resentations of the D^h point group. Then it follows that 
A ab (k ~ k') = Y^i A l ab r]i(k), and, for the frustrated ge- 



ometry of FePn, A a b(k) = A 
Cx — -cos/c^ , etc. 

Interestingly, N z ~ (r 
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) With (N a ) = 

(12k ^ab ( k ) c t,k,a c a,k,<r) ( witn a = XZ 1 V z )i which is now 



finite, is exactly associated with orbital nematicity. And 
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the degeneracy of the xz, yz bands in the T-structure, as 
seen by considering that: (i) its form factor, A^ b (k), dif- 
fers from those for the xz,yz bands (ii) in the real FePn 
structure, details of p — d orbital overlaps along a and 
b axes leads to important differences in the xz,yz band 
dispersion ll| already at LDA and Hartree-Fock (HF) 
levels. With an anisotropic form factor in (i) above, we 
get (N e ) > [13 ], as in bilayer ruthenates, already at 
HF level. Note that this is now an ON instability arising 
from an interplay between residual interactions and the 
renormalized LDA +DMFT band structure, and arises di- 
rectly from the incoherent metal. Our mechanism is thus 
radically different from one where ordered states arise 
from a weakly correlated LFL metal A finite A^$(k) 
will split the xz, yz degeneracy above and modify the 
xz, yz band dispersion anisotropically, with an increase 
in xz orbital character in the "reconstructed" FS across 
the T-0 transition, as noted recently [ll|, an effect we 
will consider in future work. The symmetry dictated cou- 
pling to the sizable magneto-elastic interaction in FePn 
will reinforce this splitting by coupling (N) linearly to 
the strain. This acts like a sizable "external" field, po- 
tentially enhances extant nematic order and will smear 
the nematic transition, making it impossible to separate 
nematicity from the T-0 distortion. 

If the xz band is thereby pushed below its position 
in the T-phase, it is not hard to deduce (see results 



below) that n x 



> n. 



giving a finite value for (N). 



Strong local correlations (U,U',Jh) will now have a 



stronger ("Mott") localizing effect on the (more popu- 
lated) xz band states, and we propose that the short- 
ening of the lattice constant b is thus a direct fall-out 
of this correlation-induced instability to the ON phase. 
This will immediately induce the T-0 distortion (b < a) 
and create conditions propitious for a q = (tt, 0) AF in- 
stability. Enhancement of Mottness in the xz sector, 
together with the incoherent "normal" state above T c , 
is also consistent with the observation of an insulator- 
like Pb(T) below 200 K in the 122-systems, and cannot 
be found in weak-coupling pictures. Thus, the observed 
ON-phase (competing with U-SC) is intimately linked to 
enhancement of the OS incoherence, already a salient fea- 
ture of the "normal" state. Sizable multi-orbital (MO) 
correlations and OS Mottness thus turn out to be neces- 
sary for even a qualitative reconciliation of these unusual 
features, as we detail below. 

Extending earlier LDA+DMFT studies, we focus on 
a four band model, since the ON phase involves pre- 
dominantly planar bands. Defining 5 — as the 
orthorhombicity, the orbital-lattice coupling is H e ; = 
^5{n xz — n y z) + , a la the more familiar Jahn- Teller 
coupling in MO systems. This self-consistently renormal- 
izes both, (N) and 6, as alluded to before. Since ab initio 
values of A, C s o are not known, we adopt the strategy of 
using trial values for (N') — (n xz — n yz )o — — in ad- 
dition to the (N e ) computed fro m a (nh.a.<j —nk,b,o-)) 
computed within LDA+DMFT |l2j. We input trial val- 
ues of (N) = (N e ) + (N 1 ) as "external" symmetry break- 
ing fields in the xz,yz orbital sector, and recompute the 
LDA -(-DMFT spectral functions in the sizably correlated 
multi-band Hubbard model [l3j, but now including the 
finite (N) above, along with U = 4.0 eV, U' = 2.6 eV and 
Jh = 0.7 eV, as used in earlier work. (N z ) is now self- 
consistently computed within DMFT from the converged 
G X z{k,w),G yz (k,uj). 

We now discuss our results. In Fig. [lj we show the 
(converged) LDA+DMFT spectral functions for the in- 
plane d orbitals with (N) = and (N) ^ 0. Clear low- 
energy OS incoherence is seen in the "normal" metallic 
(T) phase, a direct consequence of sizable MO correla- 
tions, as noticed by many workers [HI]. With (N) ^ 0, 
closer examination reveals interesting changes germane 
to the above discussion: (i) (N) > and "zeeman field" 
(h z > 0) remove the xz — yz degeneracy of the T-phase. 
With U = 4.0 eV, U' = 2.6 eV, MO Hartree shifts from 
sizable MO correlations renormalize h z = 5. DMFT val- 
ues for (N) are found to be 0.049 for x = 0.2(<S = 0.4) eV 
and 0.0315 for x = 0A(S = 0.24) eV, vanishing in the 
T phase. Concomitantly, 5(T) follows the ON order pa- 
rameter via S = —fj-^(N). The dynamical self-energies 
(S a (w)) from DMFT then drive sizable changes in spec- 
tral weight transfer (SWT) over energies O(1.0 — 2.0) eV 
in response to the renormalization of S. This feature is 
generic to strongly correlated quantum matter: obser- 
vation of sizable SWT across the T-0 transition should 



3 




m(eV) <0(eV) a>(eV) m(eV) 



FIG. 1: (Color online) Orbital-resolved LDA+DMFT (with 
U = 4.0 eV, U' = 2.6 eV and J H = 0.7 eV) density-of-states 
(DOS) for the Fe d orbitals for three 5 values (top) and three 
doping values with fixed S = 0.4 (bottom). Large-scale and 
anisotropic dynamical spectral weight transfer accompanying 
the incoherent metal-nematic metal transition is clear, in good 
accord with optics 



FIG. 2: (Color online) Imaginary parts of the orbital-resolved 
LDA+DMFT self-energies for S = 0.4, showing sub-linear 
in-u dependence for to close to Ef, reflecting an incoherent 
metal. Also, the sharp pole-like feature in lmE xz progres- 
sively sharpens and moves closer to Ef with S, reflecting en- 
hanced OS incoherence. 



put our approach on firmer ground, (ii) revealingly, the 
pseudogap (PG), already a feature of the "normal" state, 
deepens for the xz band and reduces for the yz band. 
Concomitantly, the pole-like structure in ImS IZ (aj) pro- 
gressively sharpens up and moves closer to Ep with in- 
creasing x (Fig. [2]) for not-too-large x. Finally, given 
sizable U', interband "dynamical proximity" effects in- 
duce sizable spectral changes in the xy, x 2 — y 2 spectra 
as well. Thus, OS incoherence, already present in the 
"normal" state LJ, [l4| , is further enhanced across the 
ON transition (iii) Using n a = (— 1/tt) ff^lmG a (ui)duj, 
we find n xz > n yz , as advertized above. Thus, (N) > 
and, given the enhanced (Mott-like) localization in the 
xz sector, the lattice constant b now shrinks relative to 
its value (b = a) in the T-phase: this is precisely the 
T-0 distortion. With finite (N) ^ tied to enhanced 
Mott-like localization in the ccz-sector, one can now write 
an effective ferro-nematic pseudospin model [6|: H e ff ~ 
-3 £<ij> NfN? - 8(T) J2i N?, and, with 6(T) as cho- 
sen above, structural anisotropy disappears at Tt~o- But 
strong precursor nematic correlations survive appreciably 
above Tt~o (cf. Ising model in a zeeman field). 

//incipient nematicity exists above Tt-o, one expects 
electronic anisotropy to persist up to a sizable T > Tt-o- 
this should manifest as a in-plane resistivity anisotropy 
above Tt-o- In addition to the DMFT contributions, 
we also accounted for the influence of the (short-range) 
non-local ON (plus T-O) fluctuations on the dc resis- 
tivity (these are necessary close to Tt-o) by including 



coupling of carriers to the static, T-dependent nematic 
susceptibility 15[ using an embedded cluster-variation 
method for H e ff above. We used the DMFT effec- 
tive masses and J\ , J2 as above in the evaluation of the 
short-range order contribution to p a ,b(T). This nicely 
fits within our DMFT-based view, since non-local (Ising) 
correlation effects are only included at a static level. 
In Fig. |31 we show the total dc resistivities p a (T) and 
Pb(T) calculated within LDA+DMFT-plus static non- 
local ON correlations. Quite remarkably, we find that 
the in-plane resistivity anisotropy indeed survives up to 
T ~ 250 - 300 K, sizably above T s = 150 K, though it is 
a bit overestimated at high T. It is instructive to notice 
that this anisotropy already exists within DMFT (inset 
to Fig. [3]), and is enhanced sizably near TV_o due to 
non-local nematic correlations (which peak at T^-o i n 
our H e ff), in agreement with data. Here, Pb(T) > p a (T) 
in spite of b < a (S(T) > 0) naturally follows from in- 
creased tendency to (Mott) localization in the xz band 
sector (corresponding to orbitals pointing along b) across 
the ON- plus T-0 transition (see Fig. Q]). Very satisfy- 
ingly, this anisotropy increases with indeed ob- 

served in experiment (in a 122-FePn recently, but this 
also could be seen in 1111-FePn systems, where Tt-o 
and Tsdw are separated). We emphasize that this is in 
conflict with pure J\ — J2 spin-only models, where one 
expects Pb(T) < p a (T). It is also inconsistent with a 
weakly correlated Landau Fermi liquid (LFL) view for 
same reasons. In contrast, a ON-plus T-0 transition 
within our OS-Mottness view fulfils all above experimen- 
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FIG. 3: (Color online) Resistivity versus temperature, show- 
ing the in-plane resistivity anisotropy with (main panel) and 
without (inset) the non-local, static (Ising) nematic correla- 
tions. Notice that DMFT already yields the correct sign of 
the anisotropy, and non-local corrections further enhance it 
near Tt-o, in good accord with data. 



tal constraints [2j. Finally, the optical conductivity in 
DMFT (solely involving the full DMFT Green functions, 
G xz {ui),G yz {ui)), immediately implies different optical 
conductivities along a, b. Sizable SWT, O(2.0) eV, across 
this ON-plus T-0 transition in DMFT DOS above should 
show up as SWT in o a ^(ui) across Tt~o- Remarkably, 
both these features are indeed seen [16j |. lending strong 
support for our view. 

A finite (N) > also implies ferro-orbital (FO) order- 
ing, and, with sizable magnetoelastic coupling, leads to 
conditions propituous for the Q = (tt,0) striped SDW 
order to occur. Microscopically, the "localized" orbital 
pseudospins directly follow from enhanced tendency to 
Mott localization of the xz orbital states in our work. 
The "super-exchange" couplings in the effective Ji — J 2 
spin model, with "localized" spins arising from the inco- 
herent Hubbard bands are now explicit functions of 
the neighboring (Ising) orbital correlations along a and 
b @, and, for realistic Ji a , Ju, J2, the striped AF-SDW 
state is indeed the ground state. Thus, AF-SDW order 
will always follow, or be coincident with, the T-0 dis- 
tortion, as is well known. Interestingly, the anisotropy 
of the spin-wave excitations [l^ in q-space also follows 
from such a "localized", frustrated model [3]. A weak 
coupling LFL approach Q can also achieve this upon 
"tweaking" the ellipticity of the electron pockets phe- 
nomenologically. However, it is clearly inconsistent with 
Pb{T) > p a {T), as well as with the insulator- like pb{T) 
above Tt-o- Finally, our work is not inconsistent with 
LDA+DMFT works for the AF-SDW state [H where co- 
existing quasi-itinerant and localized band states are in- 



voked; however, it bares the intimate connection between 
the ON-plus T-0 transition and subsequent AF-SDW in 
FePn. 

Thus, our work offers a unified rationalization of 
various unusual features seen recently in underdoped 
122-FePn. The observed in-plane resistive and optical 
anisotropics are understood as consequences of an elec- 
tronic ON-plus orthorhombic order arising as a p-h in- 
stability (hence, competing with U-SC [H) of the "high" 
T OS-incoherent metal. Finding of nanoscale nematic 
domains in 122 systems, as well as nanoscale electronic 
order in the 1111 systems is not inconsistent with effects 
of (doping induced) static disorder on our ON-plus or- 
thorhombic state, which falls into the "liquid gas" univer- 
sality class of an Ising model in a (random in presence of 
disorder) zeeman field [2(| • To the extent that such fluc- 
tuating ON order will induce a spin gap and two distinct 
local "environments" (corresponding to (N) = +n, —n), 
this should "rationalize" NMR data Q as well. Our work 
opens up an intriguing possibility: an ON order with a 
T-0 distortion, with subsequent stripe- AF, might be the 
primary competitor to unconventional SC in FePn, and 
establishes the important role of sizable MO correlations 
and OS Mott physics in this context. 
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